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Abstract. Using simple scaling arguments and two-dimensional numerical simulations of a granular gas 
excited by vibrating one of the container boundaries, we study a double limit of small 1 — r and large L, 
where r is the restitution coefficient and L the size of the container. We show that if the particle density 
no and (1 — r 2 )(noLd) where d is the particle diameter, are kept constant and small enough, the granular 
temperature, i.e. the mean value of the kinetic energy per particle, (E)/N, tends to a constant whereas 
the mean dissipated power per particle, {D)/N, decreases like 1/yN when N increases, provided that 
(1 — r 2 )(noLd) 2 < 1. The relative fluctuations of E, D and the power injected by the moving boundary, /, 
have simple properties in that regime. In addition, the granular temperature can be determined from the 
fluctuations of the power I(t) injected by the moving boundary. 



PACS. PACS-45.70.-n Granular systems 
noise, and Brownian motion 



PACS-05.40.-a Fluctuation phenomena, random processes, 



1 Introduction 

Granular gases have been widely studied in the recent 
years as a canonical example of dissipative systems. En- 
ergy is lost at each collision between macroscopic particles 
but a statistically stationary regime can be reached with 
a continuous energy input, achieved in most experiments 
by vibrating the container boundary. In that regime, the 
power / injected into the gas by the moving boundary, and 
the one dissipated by inelastic collisions, D, have equal 
time averages (/) = (D) and the granular temperature, 
i.e. the mean kinetic energy per particle, (E)/N is con- 
stant in time. Several laboratory experiments p] and nu- 
merical simulations [2] have been performed to study the 
properties of the granular temperature. It has been shown 
that / and D display strong fluctuations [3] and this has 
been confirmed experimentally [I] . Besides the obvious re- 
lation, (/) = (D) , there also exist relations between higher 
moments of the fluctuations of / and D and their large de- 
viations [BJ. 

Like in most dissipative systems [3], the energy bal- 
ance takes the form E = I(t) — D(t). Our aim here is to 
study a small dissipation limit such that (/) /N and (D) /N 
vanish whereas (E) /N stays finite. We show that this can 
be achieved if the restitution coefficient, r, is scaled ap- 
propriately, (1 — r 2 )(noLd) — constant when the size L 
of the container is increased at constant particle density 
no- In that quasi-elastic limit, the granular temperature 
(E)/N becomes constant whereas the mean injected or 
dissipated power per particle scales like l/^/N, provided 



that (1— r 2 )(noLd) 2 < 1. The size of the system cannot be- 
come arbitrarily large, otherwise inhomogeneities appear 
as discussed in section 2. The probability density func- 
tion (PDF) of the particle velocity becomes close to Gaus- 
sian and relative fluctuations of kinetic energy or injected 
power behave according to the law of large numbers. In 
addition, the temperature extracted from the fluctuations 
of injected power turns to be quite close to the granular 
temperature. 



2 Scaling arguments and the quasi-elastic 
limit 

A two-dimensional granular gas is simulated with an event 
driven molecular dynamics method [5\ . N disks of mass m 
and diameter d are enclosed in a square box of size L. En- 
ergy input is provided by one vibrating wall with a sym- 
metric saw-tooth motion, y = ±Voi, of period T. Like in 
most previous studies, particle- wall collisions are elastic 
whereas inelastic binary collisions between particles are 
modeled with a constant restitution coefficient r, where r 
is the ratio between pre- and post-collisional normal rela- 
tive velocities. 

When r is sufficiently close to one and the gas is suffi- 
ciently dilute, the medium is roughly homogeneous, n(y) » 
uq. Then the mean dissipated and injected powers can 
be easily estimated. (D) can be related to the mean ki- 
netic energy (E) by observing that an amount of energy 
(1 — r 2 )(E)/N is on average lost at each binary collision. 
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This suggests (D) cx (1 — r 2 )(E)v c , where v c is the col- 
lision frequency per particle, v c cx nQdv rrns whith (E) = 
mNv 2 ms /2. The mean injected power (J) has been eval- 
uated in several papers, with or without gravity 0EU2]- 
For a symmetrical excitation, it has been shown that (I) cx 
rriVftVp, where vp is the collision frequency of the particles 
with the moving boundary [8]. In the present conditions, 
vp cx noLv rms and we obtain (I) cx raV^ n§Lv rms . Then, 
using (I) = (D), we find for the granular temperature 



(E) 
N 



oc 



{l-r 2 )n Ld 



while the mean energy flux per particle is given by 
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We have shown that for r = 0.95 and r = 0.9, the above 
relations give correct predictions in the dilute limit [7]. 
Notice that the elastic limit r — * 1 at fixed L is singular 
since the granular temperature E/N diverges. In order to 
obtain a finite limit, we must take the double limit r — ► 1 
and N large with the constraint 

(1 — r 2 )(noLd) = constant = C with no = constant. (3) 

Then , the granular temperature is kept constant while the 
mean energy flux per particle from injection to dissipation 
decreases like 1/Vn. 

Those simple behaviors no longer hold if the particle 
number and the size of the container become too large, 
the density being kept constant. Indeed, for L 3> Iq, a 
particle undergoes (L/Iq) 2 collisions when traveling on a 
distance L. Thus it dissipates an amount of energy of or- 
der (n Ld) 2 (l - r 2 )(E) /N. The validity of the above scal- 
ing laws then require (1 — r 2 )(noLd) 2 < 1, otherwise the 
granular temperature and density profiles are non homo- 
geneous. For a given value of C, this provides an upper 
bound for L and N. 



3 Numerical results 

We take m and T respectively as units of mass and time, 
Vq = 3, and d = 0.5. In most simulations, tiq = 0.01 and 
C = 0.025. Then, we vary the particle number in the range 
25 < N < 2500. We observe in Fig. la that the granular 
temperature is constant within 3% in that range. Its nu- 
merical value (320) is in agreement with our order of mag- 
nitude estimate (360) from {J). In addition, as shown in 
Fig. lb, the energy flux per particle decreases according to 
(I)/N cx 1/VN as expected. Thus, the above scaling laws 
can be observed on three decades in N in the quasi-elastic 
limit ^ whereas they are restricted to less than a decade 
when N is changed with no and r kept constant [7 . This 
results from the homogeneity of the system which is main- 
tained for larger values of N in the quasi-elastic limit ([3|) . 
In addition, the PDF of the velocity fluctuations remains 
almost unchanged and close to Gaussian in that limit (see 
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Fig. 1. a) Granular temperature (E)/N (o) and temperature 
Tf — 1//3 (*) extracted from the fluctuations of the injected 
power 7 as a function of N for fixed density, no = N/L 2 = 0.01 
and C — 0.025. The inset shows that (E) /N does not reach a 
plateau if no = 0.04 and C = 0.45 b) Mean injected power per 
particle, {I)/N as a function of N, for fixed density, no — 0.01 
and C = 0.025. The full line has slope -1/2. 



Fig. 2). The velocity V along the vibration axis is less 
Gaussian than the one of the perpendicular component 
U when iV is small. This slight anisotropy is not surpris- 
ing and does not persist when N increases. Although it is 
known that the PDF cannot be exactly Gaussian in gen- 
eral [TU] , it is very close to a Gaussian in the limit we are 
considering. That nearly equilibrium behavior is also dis- 
played by relative fluctuations of energy E and power / or 
D. Their standard deviation cr(X) related to their mean 
value are shown in Fig. 3. We observe that they are in 
agreement with the predictions of the law of large num- 
bers, like for systems at equilibrium. Since E and D are 
quantities averaged on the whole surface L 2 , their relative 
fluctuations decrease like L _1 i.e. 1/VN, whereas for / av- 
eraged on the moving boundary, relative fluctuations de- 
cay like A/" 1 / 4 . Thus, we have a(E) cx y/N, a (I) cx N 1 / 4 
and a(D) = 0(1). Those simple behaviors are not ob- 
served when N is increased with no and r kept constant 

It is interesting to determine how the Gaussian be- 
havior is lost far from the present limit. To wit, we take 
C = 0.1 instead of 0.025. We obtain (E)/N roughly 4 
times smaller as expected from (JXJ) and the previous results 
remain qualitatively similar although higher moments of 
the velocity PDF's display a slightly less Gaussian behav- 
ior. We then take, no = 0.04 and C = 0.44. In that case, 
(E) /N does not reach the plateau predicted by equation [1] 
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Fig. 2. PDF's of the velocity component V parallel to the 
vibration axis, for TV = 25 (dotted line), TV = 40 (dash-dotted 
line), TV = 500 (dots) and TV = 2500 (dashed line), with fixed 
density, n = 0.01 and C = 0.025, (1 - r 2 )(n Ld) 2 < 0.0625 
and for TV = 2500 (*) with n = 0.04 and C = 0.45 thus 
(1 — r 2 )(noLd) 2 — 2.25. The full line is a Gaussian fit. 



10' 



A 10° 
V 

R 

X 10" 



10" 




10 



10 



10 



10 



N 



Fig. 3. Standard deviations of the injected power I (□), the 
dissipated power, D (*), and the kinetic energy, E (o), related 
to their respective mean values, as functions of TV", for fixed 
density, no = 0.01 and C = 0.025. The full lines have slope 
— 1/2 while the dotted line has slope —1/4. 



10" 



10" 



10" 



10" 




Fig. 4. PDF's of the injected power for m = 0.01, C = 0.025 
and different values of TV: TV = 40 (■), TV = 100 (+), TV = 500 
(o), TV = 2500 (*). The full lines are Gaussian fits. 



with the moving boundary occur during each period. Cor- 
respondingly, the standard deviation increases with TV like 
cr(I) oc TV 1 / 4 . 

We now consider the PDF, P(I T = e), of the injected 
power I(t) averaged on a time interval r for different val- 
ues of TV. In Fig. 5a, we draw centered PDF's reduced by 
the variance estimated as follows: according to large de- 
viation analysis, for r 3> tj where tj is the correlation 
time of I(t), P is of the form, P oc exp [r/(e)], where / is 
maximum for e = (I). Using the Gaussian approximation 
near (7), together with the previously observed scaling of 
the variance with TV, we find the approximate form 



P(I T 



27TCr TV 1 /4 



exp 



2al 



TV 



(4) 



as shown in the inset of Fig. la. In addition, the PDF's of 
V exhibits strongly non Gaussian behavior for the largest 
values of TV when (1 — r 2 )(noLd) 2 becomes of order 1 ((*) 
in Fig. 2). 



4 Fluctuations of the injected power 

Since the injected power I(t) can be easily measured ex- 
perimentally by recording the force applied by the par- 
ticles on the moving boundary times its velocity, it has 
been proposed to try to extract some information on the 
system by studying its fluctuations [3]. The PDF, P(I), 
of the injected power I(t) during one vibration period T 
is displayed in Fig. 4 for different values of TV. It involves 
negative events corresponding to particles colliding with 
the moving boundary at velocity —Vq. The PDF's have 
a roughly Gaussian shape contrarily to previous ones ob- 
tained outside the present scaling range, which exhibit 
asymmetric exponential tails [MTT] . This is specific to the 
limit considered here, where many uncorrelated collisions 



valid in the neighborhood of (I). In addition to scaling 
<J 2 {I T ) ex y/N/r, it is shown that the PDF's of I T for 
different values of r and TV roughly follow a single curve 
when plotted as a function of the centered and reduced 
variable (I T — (7))r 1 / 2 TV -1 / 4 (see Fig. 5a). However the 
whole distribution is skewed with a slightly larger positive 
tail, so it is not Gaussian. 

Fig. 5b shows that the quantity, - log (P(e)/P(— e)) is 
linear in e up to values of e larger than (I). Using ((4|, we 
verify that the slope is 2(1) / 'oq^/N \ which does not depend 
on t and TV (since (I) cx s/~N). The inverse of the slope is 
homogeneous to a temperature Tp = 1//3 which is plotted 
in Fig. la. We observe that T F » (E)/N within 1% in 
the range 500 < TV < 2500. This strongly suggest that the 
granular temperature can be extracted from the fluctua- 
tions of the injected power in the present limit. However, 
we emphasize that the fluctuation theorem [12] is not ex- 
pected to hold for a dissipative granular gas [3llllj. In 
general, Tp ^ (E) /TV and their respective behaviors when 
TV is varied are also different [B]. The present observation 
Tp w (E)/N is quite specific to the quasi-elastic limit ([3|). 
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Fig. 5. a) PDF of (J T - {I))/{\J 'VN '/t), for n a = 0.01, C = 
0.025 and TV = 40, r = T (+), N = 40, r = 3T (*), AT = 40, 
r = 5T (x), TV = 100, r = T (□), iV = 100, r = 3T (0), 
TV = 100, r = 5T (>), AT = 1000, r = T (o), AT = 1000, r = 3T 
(...), TV = 1000, r = 5T (•). The full line is a Gaussian fit. b) 
The quantity i log jf^^T^ , where / T is the injected power 
averaged on a time interval r, is plotted as a function of e 
(same symbols as for a). 

5 Discussion and concluding remarks 

A double limit, involving r-* 1 together with TV — * oo, 
has been first considered in [13] for a one-dimensional 
granular medium and called the quasi-elastic limit. It has 
been shown that the leading order correction to the Boltz- 
mann equation for particles without interaction is of or- 
der (1 — r)TV, so in the limit N — > oo with N(l — r) kept 
constant and small enough, inelastic collapse of particles 
is avoided. The same condition for the absence of inelas- 
tic collapse has been obtained in higher space dimension 
D, (1 — r^n^LdP" 1 small enough [T3j. The scaling used 
here (1 — r 2 )noLd D ~ 1 = constant, is only slightly differ- 
ent. However, it has been derived from the energy balance 
in an excited granular gas instead of considering the Boltz- 
mann equation for a freely cooling gas. We also mention 

that a more restrictive condition, (1— r) {n^LdP^ 1 ^ small 
enough has been given in [14] for a gas like system, i.e. non 
condensed. A condition similar to the latter can be under- 
stood in a simple way within an hydrodynamic limit. It 
can be also written (1 — r 2 ) -C K 2 where K = Iq/L is 
the Knudsen number. This is sometimes considered as a 
requirement for the validity of the hydrodynamic descrip- 
tion of granular media [T5] , although this is still a matter 
of debate [16] . 



Our simulations show that the constraint (1— r 2 )noLd D 
= constant in the dilute limit, allows us to observe a fairly 
large range of TV for which the velocity PDF are roughly 
Gaussian although the forcing remains purely determin- 
istic (no thermostat or applied fluctuating force). In that 
limit, the granular temperature remains constant as the 
size of the container is increased at constant density. The 
mean dissipated power per particle decreases like 1/s/N 
and the relative fluctuations of energy or power behave 
like in systems at equilibrium. In addition, an out of equi- 
librium temperature, defined from the fluctuations of the 
injected power, becomes equal to the granular tempera- 
ture. Those simple behaviors no longer hold if the par- 
ticle number and the size of the container become too 
large, the density being kept constant. One must also keep 
(1 — r 2 ) (noLd ^ 1 ) < 1. When that parameter becomes 
of order one, the velocity PDF is non Gaussian, the fluc- 
tuations do not scale any more according to the law of 
large numbers and the temperature extracted from the 
fluctuations of injected power becomes different from the 
granular temperature. In the quasielastic equilibrium-like 
regime, it would be interesting to determine if the non- 
equipartition of energy observed in binary granular sys- 
tems still holds [IT], or on the contrary, if other typical 
equilibrium properties are observed. 
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